load eq_SS.mat


yyaxis left
idx = find(glob.a_stat_E > 5e-7);
plot(exp(glob.sig_grid(idx)), glob.a_stat_E(idx), 'b-.', 'LineWidth', 4)
set(gca,'YTickLabel',[]);

yyaxis right
plot(exp(glob.agrid), glob.a_stat, 'LineWidth', 4)
xlim([0 6])

set(gca,'YTickLabel',[]);

legend('Signal', 'Realized productivity')

set(0,'defaultaxesfontname','cambria math') % beautify the axes a bit
set(0,'defaultTextFontName', 'cambria math')

set(gcf,'units','points','position',[10,10,500,300])
set(findall(gcf,'-property','FontSize'),'FontSize',16)


print('-dpng', 'figures/productivity_and_signal.png')